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How turbomachinery is used in Rocket Engines 



MSFC Fastrac 
engine 


• Liquid Fuel (LH2, Kerosene) and Oxidizer (L02) 
are stored in fuel tanks at a few atmospheres. 

• Turbines, driven by hot gas created by mini- 
combustors, tied with shaft to pump, which 
sucks in propellants and increases their 
pressures to several hundred atm. 

• High pressure propellants sent to Combustion 
Chamber, which ignites mixture with injectors. 

• Hot gas directed to converging/diverging nozzle 
to give flow very high velocity for thrust. 


Motivation: Avoid High Cycle Fatigue Cracking in Turbomachinery 


• Crack found during ground-test program can stop engine development 

- If crack propagates, it could liberate a piece, which at very high rotational 
speeds could be catastrophic (i.e., engine will explode) 



Structural Dynamics Basics 



Modeling and Free Vibration of SDOF Systems 


A) Equation of Motion for Undamped System 

1) Model Spring-Mass System to represent real structure 



For cantilever beam, S = 


PL 3 
3 El 


Derive k beom (in class): 

Here, P=k8 is analogous to 
typical F=kx, so k=p/S 



3EI 


2) Generate Equation of Motion (EoM) using Newton's 2 nd Law 


ku 

<=> 


m 



XF r — mil 

A 

~ ku + External = mU 


mii + ku = F 


external 





Solution to Equation of Motion (2 nd Order ODE) 


U (t) total 



particular ; steady state , nonhomogeneous 


+ 


U 


complimentary % transient , homogeneous 


Where 

• U complimentary, transient, homogeneous is the Component of the Solution when the 

RHS of the EoM is zero. 

- In physical terms, this is the response due to the internal dynamic 
characteristics of the structure, and comes about when there are 
non-zero Initial conditions (I.C/s). 


• U particular, steady state, nonhomogeneous ' S the Component of the Solution when the 

RHS of the EoM is non-zero. 

- In physical terms, this is the response due to external forcing 
functions. 


Solution Methods for Homogeneous Component 


1) simplest, worth remembering: 

— Assume solution u=u(t) is of form 


u(t) = Acos(cot) 
u(t) = -Acosin(cot) 


ii(t) = -Aar cos {cot) 

- Now plug these equalities into eq of motion: 


r\ 

m(-Aco cos cot) + k(Acosoot) = 0 
Acoscot{k - co 2 m) = 0 


For A coscot = 0, A has to = 0 , i.e., no response ("trivial solution") 
Therefore, 


Prize question: is the natural 
frequency of a system the same on 
the moon as on earth? 

Yes 

Define X = Eigenvalue = co 2— Natural Frequency 2 


2 k 
co = — 

m 



So, solution for u= u(t) 


is u(t) = Acos{ff t) 


where A depends on the 
initial conditions (1C) 



Natural Frequency Units, Period Relationship 


Question: What is relationship of the natural frequency oo (sometimes called 
"circular natural frequency) to the natural frequency/ in hz? 



Question: What is relationship of the period T to the natural frequency/ 
and the circular natural frequency oo? 


T = sec 

/ hz 


2 n 

— sec 


co 



E.G . : u = 2sin(30 

® = 3 R “ d /tc 

Q Rad/ 

f = ^ hz = 0Allhz => u = 2sin(2;r(.477)/) 

lie 

T = — sec = — sec = s = 2.094s 

/ G) 3 Rad / ec 





Damped Free Vibration of SDOF Systems 


iV 


U+ 


Damper has parameter c 



^ F external ~ 

-cu — ku = mil 
mil + cu + ku — 0 


lb 


in 

sec 


Damping Categories 


• 3 cases of solutions: 


1) Critical Damping 


c = 



critical 


2) Overdamped 


c > ijkm 

r r 
critical 

If we give a critically damped or overdamped SDOF system an 
initial displacement u 0 , we get 





3) Underdamped 

c < 2 4km 

r <C r 
^critical 

Now, to talk about 
underdamped systems, 
go back to the damped eq. 
of motion: 

mil + cu + ku = 0 
divide by m 

a + ■& u + u = o 

m m 


Now define viscous damping ratio £ as a percentage 
of critical damping 



Using co 2 = ^, we get 


” f> C cr • 

u + - c ’■ 


cr u + co u — 0 


m 


Use c cr =2 Jkm to get 


.. 2 £yjkm . 


u + 


zi + rtf w = 0 


m 


u + 2<^&> zi + of u = 0 


VERY IMPORTANT FORM, 

ALL YOU NEED ARE DYNAMIC 
CHARACTERISTICS! 




Response to Harmonic Excitation 


\ k 


N 

-1 

m 

C 



F 0 cosQt 


Q = Excitation Frequency 
f = Excitation Force (NOT enforced 
displacement or acceleration, but e.g. 
aerodynamic force) 
co = System Natural Frequency = 



eq. of motion: mil + cu + ku = F o cosQ t 


Solution: u(t) = u t (t ) + u ss (t) 

where 

u t (t ) = transient solution = homogeneous (rhs=0) sol'n 

u(t) - e~^ at ( u a cos co d t + + ^ 0M ° sin co d t) 

which decays to zero after a few cycles, 

If the excitation frequency stays constant or slowly varies. 



Steady State Solution for Non-Homogeneous Component 


Response U as a function of excitation frequency is 

U(Q) = H(Q)U st 

where we define the "Complex Frequency Response" 
Dynamic Response U 


H(Q)= 

Static Response 
and where we determine the static response U st to force F o using 


st 


kU s , = F 0 -> U s ,= y k 


After a long derivation, we get 


(l-r) +(2 Crf 


and 


i ( — 2Cr^ 

phase lag <f> = tan' 1 - r -^ y 


i-r 


j 


where we define the Frequency Ratio r = — 

co 


Resonance is defined when Q = co, ie, r=1 . 
At r=1 , | H (Q)| = ^ = Quality Factor Q 


14 



Magnification- Example 




Example: 

F=2; c=0.6; m=l; k=9 





2\[k 


m 


■= o.i 


U static — ~~T~ — 9-222 

k 


At resonance, | U | =Q U statjc = 

It/ l=— (.2222) = 1.111 

2^ 

For Q = 2.8, r= Q / = ^=.9333 

5 / CO 3 


tf(O) 


1 


^|(l-r 2 ) 2 + ( 24 'r) : 



, SO 


^ (l -0. 9 3 3 3 3 2 ) 2 + (2 *.1*0. 9 3 3 3 3)" 


= 4.408 


(j> = tan 


-2^r 
1- r 2 


= tan 


-2(.1).9333 

1-.9333 2 


_J 0.01866^1 

= tan 

^ 0.12889J 


-.9665 




Given Damped system shown below: 



W= 386.1 lb 
K= 4 Ib/in 
C= 0.3 Ib-sec/in 

C critical = 


a) What is the frequency of the steady-state dynamic 
response u(t)? 3 Rad/sgc 


e) Are there any excitation frequencies where the 
response would be close to the static response, and if so, 
what would be one of these frequencies? 


b) If the load is applied statically, what is the displacement? 

F 5 

Static =J= 4 =1 - 25 " 

c) What is the approximate maximum response of the 
mass for any frequency of excitation? 

c cr = 24km = 2^4*1 = 4 

= 0. 3/4 = 0. 075=><2 = — = = 6.667 

2C, 2*0.075 


Yes, frequencies < 0.5*co, therefore < 1 rad/s 


2) The Complex Frequency Response H is 
the ratio of the dynamic response to 

the static response 

as a function of frequency ratio r=Q/co . 



U 


dynamic 


Q u s 


= 6. 6667*1. 25" = 8. 333" 


d) At what approximate frequency does this occur? 



Multi-Degree-of-Freedom Modeling 


• Structures discretized using finite elements (rigorous mathematical 
representation of little piece of a solid) 

• Machines, other mechanical systems can be modeled fairly accurately with 
lumped parameter MDOF models (i.e., lumped rigid masses, massless springs & 
dampers). 

• Equations of Motion (EoM) 

- Newton method: easiest to use for translational systems, but very 
difficult for rotational motion 









Modal Analysis - Obtain Natural Frequencies and Modes 



Solutions for Undamped, Free Vibration of MDOF Systems with N dot’s. 

[m]{m}+[*:]{ m } = {0} 


Assume solution of form 

{u} m = {</>} m e i ^ t+ ^ 

m = 1,M where M<N (can choose to "use" less than N solutions) 
where m spatial solutions = Eigenvectors=Mode Shapes 


Discrete MDOF modes 



Alternate nomenclature: (/) ijr where / is dof, j is mode number 






Example: Natural Frequencies and Modes of Axial Vibration of Cantilever Bar 


2 " 


1 " 


960" 


E=2.9e7 lb/in 2 


A) Discretize into 2 finite elements, draw coordinate system: 



u i 


u- 


U: 




X 


m. 


L 

k i 


Write Eq's of Motion: 


m 


m 


3 ^ ul 


Lump mass at nodes: 
So 

m 1 =pAL/2 

m 2 =pAL 

m 3 =pAL/2 


m 


element 


=p A L 


ml | M 2 ) 


— —u 2 ) = m jMj 
mjMj - k x u 2 + = 0 

m { ii { + - k x u 2 - 0 


Element stiffness: 

k 1 =k 2 = AE/L 


k l (u l - u 2 ) - k 2 (u 2 - u 3 ) = m 2 U 2 
m 2 u 2 + k { u 2 — k x u x - k 2 u 3 + k 2 u 2 — 0 
m 2 u 2 — k x u { + (k { + k 2 )u 2 — k 2 u 3 — 0 


k 2 (u 2 - u 3 ) = m 3 u 3 
m 3 u 3 - k 2 (u 2 — u 3 ) = 0 
m 3 u 3 - k 2 u 2 + k 2 u 3 = 0 



Finite Element Method Derivation p. 3 


Equation of Motion in Matrix Form: 


pAL 

2 


1 

0 

0“ 


r ^ 


1 






AE 


0 

2 

0 

< 

U 0 

r + 

-1 

A 

A 

HI 


z 

L 

r\ 

u 

u 

l 


V J J 


U 




r 'n 

C) 



-1 

i 

< 

^2 > 

I 


Mo 


o 

< 0 > 


Apply Boundary Conditions (BC's): Since u 3 =0, can cross out 
corresponding row & column. 


Write system matrix [D]=([K]-oo 2 [M]), let ^=co 2 



1 

-1 

- x pAL 

“1 

0 ~ 

A 

L 

-1 

2 

2 

_0 

2 

/ 


4) Divide through by AE/L; 
Let 


M = 


pAL t2 

2 E 


1 — ju -1 

-1 2-2 ju 


then [D] = 


5) Set | D | = 0 
|1 ~M - 1 


= 0 


| -1 2 - 2 ju \ 

(l-//)(2-2//)-l = 0 


2 ju 1 -4// + l = 0 


6) Solve for roots p m \ solve for co n 

= 1 + 


4 ± Jl 6 - (4 * 2 * 1) , V2 

ju, 2 = = 1 ± — 

12 2*2 2 


//, = .2928 =^> \ — 


2 E 


T b 

pi} 0.1 


2(2. 9 el 4) 

— ^-r0.2928 = 737.1 => ax =27.15 

lb seer / /i oaj„\2 1 


=1.707 


•% = 


in 

2(2.9e7^) 


■(480i«)‘ 


, T 1.707 = 4297.5=><» 2 =65.5fi 

Wsqc i a onj^\2 z sec 


0.1^(480my 


sec 


7) Plug //, in to [d]{^} = 0=>{^}, use "max" normalization 

: (l-0.2928)^ n = 0 

.7071 -^ 21 = 0 => (f) 2X = .7071 

Let =1 =>^j = -7071 => {^} 2 = { 7 q 71 1 


2nd equation in [£>]{^} = 0 gives same thing (not linearly independent) 


/u 2 : (1-1.7071)$ 2 ~(j) 22 =0 
-.7071 (j) n = (j> 22 

Let ^ 12 = 1 => (f) 22 = -.7071 => {^} 2 = | 7 Q 7 i 


,.[o] 


1 1 
.7071 -.7071 


Max normalized mode shapes are useful for visualization 



8) If desired, obtain Mass Normalized Mode Shapes 
s x = 0.10206 and s 2 =0. 10206, so 



mass 


w Mode 1, both dof's IN-PHASE 

— Mode 2, both dof's 
— OUT-OF-PHASE (slinky) 

".10206 .10206 " 

.07216 -.07216 


Natural Frequencies, Modes, & Modal Matrix 


See great animations of MDOF systems by Dr. Dan Russell, Graduate Program in 
Acoustics, Penn State. Or https://www.youtube.com/watch?v=kvG70riBirl 

Eigenvalue X = natural frequency co 2 


^2 — {^}2 
/I3 = co 3 => |^} 3 



Displacements for all 
locations of mode 
shape are either in- 
phase or 180° out-of- 
phase with each 
other, but have phase 
relationship of otj with 
excitation. 


Modal Matrix <D = 


M, W 2 H! 


0 n $3 

</>nl • • • ^n3 


23 



Application of Structural 
Dynamics to Turbomachinery 


Characterization of Excitations - Speed Range 




First obtain speed range of operation from 
performance group. 

- For Rocket Engines, there are generally several 
"nominal" operating speeds dependent upon 
phase of mission (e.g., reduce thrust during "Max 
Q"). 

- However, since flow is the controlling parameter, 
actual rotational speeds are uncertain (especially 
during design phase) 


Rated Power 



Level 

70% 

100% 

Low Range 

20759.4 

26125 

Nominal 

21852 

27500 

High Range 

22944.6 

28875 


— For new LPS engine being built at MSFC, assuming 
possible variation +/-S% about each of two 
operating speeds. 

In addition, speed generally 
isn't constant, but instead 
"dithers". + 

implementation of Speed Variation in the Structural Dynamic 
Assessment of Turbomachinery Flow Path Components 

Andrew M. Brown, R. Benjamin Davis and Michael K. DeHaye 
J. Eng. Gas Turbines Power 135(10), 102503 (Aug 30, 2013) Paper No: 

GTP-13-1206; doi: 10.1115/1.4024960 





Quantify Engine Forces using Fourier Analysis 




• Forces are not, in general, perfect sine waves (although sometimes they're close) 

• We can deal with these in two ways: 

• Represent forces as sum of Sines (Spectral, Frequency Domain Approach), sum 
response to each Sine 

• Calculate response to actual temporal (time history) loading using "impulse response 
function" 

• Spectral Approach: given a periodic but non-harmonic excitation 




Fourier Analysis 


Jean Fourier realized we can write loading p(t) as sum of average, cosines, & sines: 

00 

p(t) = a D + I [ a n cos(nQ/) + ^„ sin ( hQ/)] 

where 


n - 1 


r +7"' 

a 0 = — f p(t)dt = avg value of p(t) 

1 r 

=— f p(t) co s(nQ l t)dt 

T i 

b n = ~\ P(0s ininClfidt 

Textbook Example: Using Fourier Series, represent square wave excitation as: 



p(t) = —^~ ^ f— sin^QjO 

K n=l,3,5,... \ n J 




Characterization of Mechanical Excitation due to Unbalance 



• Since all real machinery has some level of unbalance and shaft whirl, sinusoidal 
("harmonic") loads at 1 and 2 times rotational speed ("1 and 2N") will be 
generated, along with up to their 3 rd multiples (also called "harmonics") ^ 
3-6N throughout turbopump. 

Space Shuttle Main Engine Powerhead Cross-Section 



Chara 



Nozzle (37) 

Convergent-Divergent Channel 
Machined Vanes 


Manifold 

Radial inlet Jonas 


Cast Blades 


Two-Row Velocity 
Compounded Turbine 


1 st & 2n/stage Disks ° r (57J 

Cuntf Coupling uasl vanas 


Exit Guide Vanes (52) 

Zero Swirl at Discharge 


r , 


1st Row Blades (69) 
2nd Row Blades (59) 

Zero Reaction 


of Fluid Excitation 


• Harmonic excitation at engine order = 
Number of flow distortions and up to 
their 3 rd multiples arising from adjacent 
upstream and downstream blade and 
vane counts. 


• Use CFD to generate Loading 






Engine Example of Application of Fourier Series 


• Excitation wave based upon a pump with 3 primary distortions (e.g. diffusers), 
within slightly asymmetric overall field. 


• Let primary excitation at 3N have an amplitude of 1, and asymmetric primary 
distortion have an amplitude of 0.1. 


— each of these will have a harmonic, since they aren't perfect sinusoidal distortions, 
such that the harmonic of the asymmetric is 0.05, and the amplitude of the 
harmonic of the primary distortion is 0.25. 


- So p(t) = b x Sin(i) + b 2 Sin(2t ) + b 3 Sin(3t ) + b 6 Sin(6t ) 


Where 

b^O.l, b 2 =0.05, 
b 3 =1.0, b 6 =0.25 

• Have to assess 
dynamics for each 
frequency component 
of excitation. 



Now Structure: Create FEM of component. Modal Analysis 


Example: 

Turbine 

Blades 





Mode 12 at 36850 hz 


Mode 13 at 38519 hz 



ANSYS 11.QSP1 
JUL 31 2009 
11:34:15 
j2f Ibid 28k 
PLOT NO. 139 
NODAL SOLUTICN 
STEP=2 
SUB =1 


TIMB=2 

USUM 

RSYS=0 

DMX = . 247E-03 


Sm = . 254E-05 
SMX =.247E-03 


□ 

□zi 


. 254E-05 
.297E-04 
. 569E-04 
.841E-04 
. 11 IE-03 
•138E-03 
. 166E-03 
. 193E-03 
. 220E-03 
.247E-03 


dia #=5 

spd=29879, (5 of 10) 
freq ( force) =36850 
Fourier coeff=74 



ANSYS 11.0SP1 
JUL 31 2009 
11 ; 
j2f_lbld_28k 
PLOT NO. 149 
NODAL SOLUTION 
STEP=6 
SUB =1 
TIME=6 
USUM 
RSYS=0 

DMX =.160E-03 
.366E-05 
. 160E-03 
. 366E-05 
.210E-04 
.383E-04 
. 557E-04 
.730E-04 
. 904E-04 
.108E-03 
. 125E-03 
.142E-03 
.160E-03 


SNN = 
SMX ■ 


□ 

[=3 


i i 


flats 


dia #=5 

spd=31232 / (6 of 10) 
freq (force) =38519 
Fourier coeff=74 


Modal Animations 
very useful for 
identifying 
problem modes, 
optimal damper 
locations 






hi 


Create “Campbell Diagram” 


• Simplest Version of Campbell Diagram is just a glorified Resonance Chart. 


J2-S Fuel Turbopump 1st Stage Turbine Blade Campbell Diagram 


10 < ND=5; Excitation Order = 74 (2x37 Nozzles) 



2 8 1 1 1 1 1 1 1 1 1 

' 2.4 2.6 2.8 3 3.2 3.4 3.6 


Pump Speed (RPM) 



Modal Analysis has Multiple Uses 


• Redesign Configuration to move excitations ranges away from 
natural frequencies 

• Redesign component to move resonances out of operating 
range. 

• Put in enough damping to significantly reduce response 

• Use as first step in "Forced Response Analysis" (applying forces 
and calculating structural response). 


Thickness (in) 



Patran 2012.2 64-Bit 24-Apr-14 07:59:07 3.48+001 

Fringe: JULY 8 MOD SECTOR STATOR DESIGN ND30, Mode 3:Freq. =38969., Eigenvectors, Tran slatg)ggl + ^^| 
Deform: JULY 8 MOD SECTOR STATOR DESIGN ND30, Mode 3:Freq.=38969., Eigenvectors, Translational, 


ir-h/DC i/ n T'V/W 


1.86+001 


1.63+001 


1.39+001 


1.16+001 


9.29+000 


6.97+000 


4.65+000 


2.33+OOOHI 
8.50-003® 
default_Fringe : 

Max 3.48+001 @Nd 1438 
Min 8.50-003 @Nd 2766 
default_Deformation : I 


Max 3.48+001 @Nd 14381 


LPSP Turbine Stator Redesign to Avoid Resonance 


Modal analysis of original design indicated 
resonance with primary mode by primary 
forcing function. 

_ Since excitation simultaneously from 
upstream and downstream blades, critical 
to change design to avoid resonance. 

_ Extensive optimization effort performed to 
either move natural frequency out of 
range and/or change count of turbine 
blades to move excitation. 


Stator Airfoil Thickness Changes 


R03b2 


R03b2 t2 


R03b2 t4 


R03b2 t5 


R03b2 t6 


% Airfoil Span 










gig 


5 


I 

5k 










Final and Original Campbell of Modes for Stator 
Vane 30ND .Family 



Range of +/- 
5% on 
natural 
frequencies 
to account 
for modeling 
uncertainty 


LPSP Lox Pump Turbine Stator Campbell 








Damping 


A. Brown 
MSFC Propulsion 
Structural Dynamics 
36 

• Damping is critical parameter for forced response prediction, so 
"whirligig" test program used to obtain data. 

• Whirligig is mechanically-driven rotor with bladed-disk excited by 
pressurized orifice plate simulate blade excitation. 

• Key assumption is that this reflects true configuration. 

• SDOF Curve fit 
technique applied to 
selected top- 
responding blades 
to derive damping 
from response. 





Damping Results from Whirligig 


Data shows wide-variation in 
damping, but reasonable 
population (15-20 acceptable 
samples) for statistical 
characterization. 

Lognormal distribution fits 
obtained for each mode. 


Amp 


Freq 


Zeta 


Nodal Diameter 

Mode 

Samples 

Mean 

Sigma 

Min 

Max 

Mean 

Sigma 

Min 

Max 

Mean 

Sigma 

Min 

Max 

LogNormal Dist.: 

0a Equivalent 
-a Equivalent 
-2a Equivalent 
-3a Equivalent 


5 

3 

18 

15.6 

3.2 

9.9 

20.3 

10967 

17 

10936 

10997 

0.404 

0.103 

0.314 

0.720 


0.391 

0.305 

0.237 

0.184 


5 

4 

17 

7.8 

1.9 

5.0 

11.2 

13831 

69 

13695 

13908 

0.702 

0.163 

0.520 

0.976 


0.684 

0.544 

0.433 

0.343 


23068 

282 

22921 

23816 

0.146 

0.023 

0.106 

0.191 


0.144 

0.123 

0.105 

0.090 


0.183 

0.132 

0.095 

0.068 


0.223 

0.149 

0.099 

0.066 


A, Brewn 
MiFC Prepuliien 
itruetunl Dymmiei 




Can Also Use Modal Analysis in Failure Investigations 



• Examination of Modal Stress Plots provides link to location of observed cracking. 






Now, if resonance, forced response required, need to know about 

Generalized Coordinates/Modal Superposition 


• Frequency and Transient Response Analysis uses Concept of Modal 
Superposition using Generalized (or Principal Coordinates). 

[M]{u} + [C]{u} + [K]{u} = {P(t)} 

• Mode Superposition Method - transforms to set of uncoupled, SDOF equations 
that we can solve using SDOF methods. 

• First obtain [0] mass . Then introduce coordinate transformation: 

M 

| w } = n[0]{ti}m 

[M][O]{^} + [C][<I>]{i7} + [^][®]{7} = {P(0} 

I 

[^]W + E^]W + ['*\]W=[°] r { i> ( r )}- 



for the SDOF equation of motion, 

• • • 2 

mii + cu + ku = F — » ii + + 6) u = F 



So we get the same equations in q: 

iim + 2 £» m + 4 a = {</>}[ {^(0} 



• For "Frequency Response" Analysis, 
apply Fourier coefficients coming from 
CFD such that excitation frequencies 
match Campbell crossovers. 


Maximum Principal Stress 



Q 

co 


SSME HPFTP 1 st Blade Frequency Response 


6.00E+04 


5.00E+04 


4.00E+04 


3.00E+04 


2.00E+04 


1 OOE+04 


O.OOE+QO 

30,000 35,000 40,000 45,000 50,000 55,000 €0,000 65,000 



Frequency (Hz) 




Forced Response Analysis in Failure Investigations 




Cyclic Symmetry in Turbine Components 


• Many structures possess some kind of symmetry that can be used to 
simplify their analysis. 

• A cyclically symmetric structure possesses rotational symmetry, i.e., the 
original configuration is obtained after the structure is rotated about the 
axis of symmetry by a given angle. 

• Instead of modeling entire structure, only model one sector. 

• For turbomachinery structures, structural analysis is generally only possible 
by taking advantage of huge reduction in model size by using cyclic 
symmetry. 




Section courtesy of 
Dr. Eric Christensen , 
DCI Inc. 





• Most Nodal Diameter modes exist in pairs, 
same shape but rotated by ti/ND 

• For purely cyclic symmetric sections, 
highest Nodal Diameter possible is N/2 for 
even # sectors, (N-l)/2 for odd # sectors. 

• First family of modes: 

- Has unique eigenvalues 

- Has unique eigenvectors 

- All segments have same mode shape 

• The next family of modes: 

- Pairs of degenerate (identical) 
eigenvalues 

- Non-unique eigenvectors 

• The last family of modes: 

- Only exist if N is even 

- Has unique eigenvalues & eigenvectors. 



Flat Plate 
Example 
N = 8 
Segments 






Example - SSME Fuel Turbopump 3 rd Stage Impeller 


Cyclically symmetric sections have N = 6, so max ND=N/2=3 for 


those sections. 

However, much of impeller is disk 
with axisymmetry, having infinite 
# ND's. 




• Codes generate "Harmonic Families" 
of modes which only contain Nodal 
Diameter modes according to 
following stair-step pattern: 

• HO contains ND= 0, 7, 14, 

• HI contains ND= 1, 6 ,8, 13, 15,... 


• H2 contains ND= 2, 5, 9, 12, 16 


• H3 contains ND= 3, 4 ,10,11, 17 


Mode 

Natural Frequencies (Hz) 

Harmonic 0 

Harmonic 1 

Harmonic 2 

Harmonic 3 

1 

3,710 

4,210 

5,180 

4,900 

2 

8,360 

4,210 

5,1 80 

6,530 

3 

9,780 

7,190 

5,750 

6,850 

4 

10,100 

7,1 90 

5,750 

10,400 

5 

10,500 

8,740 

7,900 

11,500 

6 

11,200 

8,740 

7,900 

13,200 

7 

12,100 

9,670 

10,000 

14,300 

3 

13,500 

9,670 

10,000 

14,900 

9 

16,200 

11,200 

11,500 

15,200 

10 

17,500 

11,200 

11,500 

15,900 

11 

17,800 

12,500 

13,300 

16,000 

12 

1 8,000 

12,500 

13,300 

16,500 

13 

18,600 

13,000 

14,000 

18,900 

14 

1 8,900 

13,000 

14,000 

19,300 

15 

19,500 

14,700 

14,400 

19,600 


• • • 





Implications of Cyclic Symmetry - Generalized Force 


Generalized (or Modal) Force defined as 

{ fL= ff>. 

This is just the dot product of each mode with the excitation force vector 
and means that the response is directly proportional to the similarity of the 
spatial shape of each mode with the spatial shape of the force. 

For pure harmonic waves, the "Orthogonality Principle" states 

r f ^whenn=m 

[ sm(nO)sm(mO) = { 

1 0 otherwise 

-71 

Think of the {O} as a continuous function, and the force the same way. 

_ Then the Dot Product is the same as an integration of the product of the two 
functions. 

_ So this says that the only non-zero result of an excitation wave shaped like a 
Sine and a mode shaped like a Sine is for the components of those waves that 
have the same wave number! 


Have to determine Nodal Diameter of 
Modes to identify Resonance 


• For disks and disk dominated modes, 5ND Traveling Wave will excite a 5ND mode 



5ND standing 
wave mode of 
Impeller (modal 
test using 
holography) 


5ND travelling wave Mode of Bladed-Disc 




default_Fringe : 

Max 2.81-003 @Nd 375 
Min -2.41-003 @Nd 1011 
default_Deformation : 


Max 2.81-003 @Nd 375 
Frame' 1 


• On the other hand, 3ND 

excitation (perhaps from pump 
diffusers) will not excite a 5ND 
structural mode. 









Impeller Campbell Diagram 


• "Triple Crossover 
Points" (speed, co, and 
ND) needed for 
resonance of pure 
shroud (disk) modes. 

• ND mode at exact 
spatial number of 
distortions in 
excitation. 





“Blade/Vane” Interaction causes different ND excitation 


• Sampling by discrete number of points on structure of pressure oscillation 
results in spatial Nodal Diameter excitation at the difference of the two counts. 

• E.g., a 74 wave number pressure field (coming from 2x37 vanes), exciting 69 
blades results in a Nodal Diameter mode of 69-74=-5, where sign indicates 
direction of traveling 5ND wave ( /plot courtesy_Anton Gagne) . 



Direction of 
resultant wave 
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Spatial Nodal Diameter 5 Mode Sampled by 69 blades 
from 74 Vane Pressure Oscillation 
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Theta location (radians) 



Tyler-Sofrin Blade-Vane Interaction Charts 


• Chart identifies Nodal Diameter families that can be excited 


Downstream 

74 111 148 Stator Multiples 57 114 171 

-5 N/A 
N/A 27 
N/A N/A 


• All modes in Campbell have to be from these families 

— E.g., Nodal Diameter 5, blade mode 3 (torsion) 

• Temporal Frequency of Excitation is at the engine order of the 
distortion. 

• Much of chart is marked "N/A - not applicable" because 

— Highest number of ND waves in a cyclic symmetric structure is 
N/2 or (N-l)/2 


Blade multiples 

N/A 69 12 N/A N/A N/A 

-10 138 N/A 24 -33 N/A 

N/A 207 N/A N/A N/A -21 


228 


Upstream Nozzle 
Multiples 

37 

Blade multiples 


69 

32 

138 

N/A 

207 

N/A 





Example - LPSP Turbine Blisk Aliasing Tables, 
Non-P roblematic Modal Evaluation 

• 123 blades allow (123-1)/2=61 Nodal Diameters 



1st Stage 
Nozzles 


2nd Stage Stators 

Exit Guide Vanes 


87 

174 

261 

93 

186 

279 

97 

194 

291 

blades 










123 

36 

51 

138 

30 

63 

156 

26 

71 

168 

246 

159 

72 

15 

153 

60 

33 

149 

52 

45 

369 

282 

195 

108 

276 

183 

90 

272 

175 

78 

0 

87 

174 

261 

93 

186 

279 

97 

194 

291 


• Many crossings judged low risk, or acceptable risk for non-flight program 

- 15ND, 33ND, and 45ND crossing modes eliminated due to probable low 3X forcing 
function, high frequency 

- 52ND modes extremely complicated, high frequency 

- Many modes eliminated due to non-adjacency of forcing function 

• Also have to consider mechanical excitations order 1-6N 

• In this case, no 1-6ND modes have crossings with appropriate forcing function. 




Frequency (kHz) 




Turbine Blisk Campbell for Problematic Modes in 100% PL Range 









Natural Frequencies 


ND 

mode pair 
number 

70% PL 

100% PL 

mode shape description 

excitation source and 
order 

potential solution 

36 

5 

28955.47 

28964.54 

1 st Torsion NOTE 
12/31/13- 2 nd blade -are 
nozzles problem? 

87=1 x 1st st 
upstream nozzles 

change # nozzles to 89 to 
move lower bound of 70% 
above mode. 

51 

9 

58040.23 

58048.62 

1st blade chordwise 2nd 
bending 

174= 2 x 1st st 
upstream nozzles 

change # nozzles to 89 to 
move lower bound of 70% 
above mode. 

51 

11 

75845.08 

75878.04 

1st blade spanwise TE 2 
wave 

174= 2 x 1st st 
upstream nozzles 

CFD to determine 
magnitude of 2x 

60 

12 

78335.00 

78349.16 

2nd blade bending 

186 = 2 x upstream 
2nd st stator 

CFD to determine 
magnitude of 2x 

60 

13 

79511.67 

79524.37 

2nd blade ? 

186 = 2 x upstream 
2nd st stator 

CFD to determine 
magnitude of 2x 

51 

14 

85944.41 

85972.20 

1st blade TE spanwise 1 
wave 

174= 2 x 1st st 
upstream nozzles 

CFD to determine 
magnitude of 2x 

60 

14 

86007.76 

86035.52 

1st blade TE 1.5 wave 

186 = 2 x 

downstream 2nd st 
stator 

CFD to determine 
magnitude of 2x 

60 

16 

92827.38 

92835.74 

1st blade TE 0.5 wave 

186 = 2 x 

downstream 2nd st 
stator 

change # stators to 92, 
mode will be above range 






















Spatial Fourier Analysis helpful to identify ND number 

of both excitation and modes 


• Let's say we have measured 
pressure field 


p(t,0) 

that has the unknown 
temporal and spatial form: 

5 = Sin[3t 4- 20] 4- Cos[4t 4- 30] 


■ 



• First perform Temporal Fourier 
Analysis at each value of 6, using 
Complex Form (more efficient than 
harmonic form) 

oo 

p(t)= X 


where 





Spatial Decomposition of Each Temporal 
Fourier Component 


• Now look at obtaining spatial components for each of the temporal fourier 
components ("bins"). 


Pi ) = Z ] 



where 

p(d)e ina ' e dd 



Wave number 


2-D Fourier Transform Shows Spatial Complexity of J2X turbine flow field 
and response, used in evaluation of forced response methodologies 



57 64 


Frequency Component 



Mistuning 


• Due to manufacturing tolerances, the turbine blades on a bladed disk will never 
be identical 

• Result is that the bladed-disk modes (e.g., ND17, 1 st Torsion) will have slightly 
varying natural frequencies and mode shapes; typical variation is +/- 1.5%. 

• Variation itself called Mistuning, which generates two effects in addition to 
bifurcation of individual modes into n*individual modes: 

— Localization - mode shapes warp such that maximum deflection is at a single 
location, rather than at every high point in a "tuned" nodal diameter mode. 

— Amplification - most important effect - the maximum resonant "mistuned" 
response is frequently up to twice as much as the "tuned" maximum resonant 
response. 

• Probabilistic analysis techniques required since every bladed-disk will be 
different 

• Innumerable papers and Ph.D. theses have been devoted to this topic over the 
last 40 years. 

• Tractable techniques for predicting level of amplification for a design have only 
existed since 2004. 


Localization due to Mistuning 



Ref: Rao, J.S., Mistuning of Bladed Disk Assembies to Mitigate Resonance, Altair Engineering, 2006 







Amplification due to Mistuning 


Maxiumum amplified blade is not blade with most mistuning, appears to be most 
pronounced near locations of "eigenvalue veering" on Nodal Diameter plot. 
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Can predict amplification for a given design 
using "MISER" by assuming a std deviation of 
frequency mistuning and generating a nodal 
diameter diagram for modes of interest. 



2.5 


2 - 


0.5 


Amplification 
Results for Monte 
Carlo Analysis of 
10,000 Bladed-Disks 


1.5 


£.5 


EO=74, ND=5, Case 5752 



Blade Number 


Ref: Castinier, M. P., Pierre, C., Modeling and Analysis of Mistuned Bladed Disk Vibration: Status and Emerging 
Directions, Journal of Propulsion and Power, Vol. 22, No. 2, March-April 2006. 



Conclusion 


• Structural Dynamic Analysis of Turbomachinery is critical aspect 
of design, development, test, and failure analysis of Rocket 
Engines. 

• Process of Analysis starts consists of modeling, modal analysis 
and characterization, comparison with excitation field, and 
forced response analysis if necessary. 

• Thorough understanding of Fourier Analysis, Vibration Theory, 
Finite Element Analysis critical. 

• Knowledge of Turbomachinery Design and Fluid Dynamics very 
useful. 



